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I.  INTRODUCTION 

This  report,  the  first  quarterly  technical  summary  report 
issued  by  the  SDL,  covers  the  period  April  1965  through  Juna  1965 
Previously,  reports  of  this  type  were  issued  on  a  semi-annual 
basis;  and  seven  such  semi-annual  reports  were  issued  for  the 
period  October  1961  through  Marcn  1965. 

Analyses  completed,  for  which  results  have  been  reported, 
are  discussed  in  Section  II  under  descriptive  headings.  Work 
currently  in  progress,  which  has  not  reached  the  stage  where 
results  are  available,  is  discussed  in  Section  III.  Section  IV 
contains  a  discussion  of  the  support  and  service  tasks  performed 
for  in-house  projects  and  for  other  VELA  UNIFORM  participants. 
Work  previously  completed  and  reported  is  mentioned  only  as  it 
relates  to  analyses  in  progress  during  this  reporting  period. 

Appendix  A  is  a  listing  of  those  organizations  receiving 
SDL  data  services  during  this  period;  and  Appendix  B  contains 
selected  SDL  reports  which  are  representative  of  the  types  of 
analyses  made  by  the  Seismic  Data  Laboratory. 

II.  WORK  COMPLETED 


A.  Seismic  Partial  Coherency  Study 

This  report  present  s  the  results  of  a  preliminary  inves 
tigation  into  che  application  of  partial  coherency  techniques  to 
the  problem  of  processing  seismic  data.  Ultimately,  these 
techniques  will  be  applied  to  the  detection  of  the  presence  or 
absence  of  unidentified  components  in  seismic  noise,  determina¬ 
tion  of  filter  responses  between  two  time  series,  and  description 
of  seismic  noise  fields. 

A  specific  application  of  partial  coherencies  to  seis¬ 
mology  would  be  the  modeling  of  a  seismic  noise  field.  If  a 
noise  field  has  only  one  component  propagating  across  an  area, 
then  the  output  of  one  seismometer  should  suffice  to  predict  the 
output  of  a  second.  If  there  are  two  components  in  the  field, 
then  two  element  outputs  are  required  to  predict  a  third.  If 
there  are  n  components,  n  elements  are  required  to  predict  the 
output  of  one  additional  element.  The  methods  of  partial 
coherencies  can  be  used  to  determine  when  sufficiently  many 
inputs  are  being  used  to  predict  an  output.  Thus,  it  will  be 
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possible,  by  using  partial  coherencies,  to  determine  the  number 
of  major  noise  components  present  and  the  minimum  number  of 
elements  required  to  study  the  noise  field. 

Before  the  development  of  a  general  program,  it 
was  necessary  to  study  partial  coherence  techniques  applied  to 
the  solution  of  relatively  simple  problems  under  controlled 
conditions . 

The  application  of  power  spectra  and  cross-spectra  to 
determine  frequency  response  functions  for  simple  linear  systems, 
where  a  single  input  and  output  are  clearly  defined,  is  well 
established.  Application  of  these  ideas  to  more  complex  systems 
such  as  the  earth  where  there  are  many  input  and  output  points 
is  not  so  well  known. 


Cc.' aider  the  case  of  a  simple  linear  system  with  a  well- 
defined  single  input  (excitation)  point  and  single  output 
(response)  point.  If  the  input  x(t)  is  a  stationary  random 
process  with  zero  mean,  then  the  output  y(t)  is  also  a  station¬ 
ary  random  process  with  zero  mean.  Now,  if  G  (f)  and  G  (f)  are  the 

x  y 

one-sided  power  spectral  density  functions,  and  G  (f)  the 

xy 

cross-power  spectral  density  for  the  input  and  output;  then  the 
frequency  response  function  for  the  linear  system,  H(f)  is 
completely  determined  by  the  two  relations 


G  (£)  =|H(£)| 


Gx(f) 


G  (f)  =  H(f)  G  (f) 
xy  x 


A  quantity  of  particular  importance  in  more  complicated 
situations  is  the  coherence  function  which  is  defined  by 


•5  G  (f)| 

v  2  If) - ZX — ! — 

Y  xy  lt'  Gx(f)  Gy(f) 


It  is  easily  shown  that  0  -  y  (f)  -  1,  and  for  a  linear  system , 

xy 


v  (f)  -  1 .  Thus,  the  coherence  function  may  be  thought  of 
xy 

as  a  measure  of  linear  relationship  in  the  sense  that  the 
function  attains  a  theoretical  maximum  of  unity  for  all  frequencies 
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xn  a  linear  system.  Therefore,  if  the  coherence  function  is 
less  than  unity,  one  possible  cause  might  be  the  lack  of 
complete  linear  dependence  between  the  input  and  output. 

Another  application  of  coherence  functions  to  single 
input-single  output  systems  is  to  determine  the  effect  of 
measurement  noise  on  frequency  response  function  measurements. 
The  coherence  function  in  this  case  is 


7  (f)  = 

xy 


l+fi.+e 

i  o  10 


<  1 


(4) 


where  3  .  is  the  measurement  noise  to  signal  ratio  at  the  input 
and  P  is  the  corresponding  ratio  at  the  output.  Hence,  if 
the  coherence  function  is  found  to  be  significantly  less  than 
unity,  one  possible  cause  might  be  that  the  measurement  noise 
effects  are  not  negligible  and  must  be  taken  into  account  in  de- 
termii  j.ng  the  frequency  response  function. 

Coherence  functions  are  useful  aids  in  the  analysis 
of  single  input-single  output  linear  systems.  However,  the 
area  of  major  application  is  in  the  analysis  of  multiple  input- 
multiple  output  linear  systems  which,  for  example,  could  represent 
the  response  of  the  earth  to  internal  excitations  as  measured  at 
several  different  points.  This  would  be  the  case  when  an  array  of 
seismic  instruments  is  employed.  If  it  is  assumed  that  the  earth 
has  a  linear  response,  and  that  measurement  noise  is  negligible, 
then  a  low  coherence  function  between  two  points  will  serve  to 
indicate  the  presence  of  other  factors  which  contribute  to  the 
output  but  are  not  being  considered. 

In  order  to  investigate  partial  coherence  techniques 
under  controlled  conditions,  two  demonstration  cases  were  devel¬ 
oped.  These  are  described  briefly  below 


Case  1.  Given  x^t),  *2(t),  and  y(t)  =  a^x^t)  +  a2x2(t-T). 

Determine  a..  ,  a  ,  and  T  in  the  case  where  x  ,  x 

«*•  ~  X  z 

are  independent  and  in  the  case  where  x^,  x2  are 
correlated. 
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Case  2. 


Given  x^(t) ,  x2(t),  and 

y(t)  =  a1x1(t)  +  a2x2  (t-T)  +  . 

Determine  a^,  a2?  T,  and  the  presence  of  x^ 

under  various  conditions  of  dependence  and 
independence  between  x^,  x2#  and  x^. 

Consider  the  linear  system  sketched  below 


z(t ) 


Y  (t ) 


Two  independent  seismic  sources,  x^(t)  an<*  z(t)  are  combined 
to  form  a  new  process,  x2(t),  such  that 

x2(t)  =  0.4x^(t)  +  0.6z(t)  (5) 

Secondly,  x  (t)  is  delayed  in  time  by  0.2  second  relative  to 
x^t),  and  finally  combined  with  x^(t)  to  form  y(t).  Thus 


y (t)  =  x1(t)  +  x2(  t  -  0.2)  (6) 

If  the  transfer  functions  between  the  two  inputs  and  the  output, 

denoted  by  H1w(f)  and  H_w(f) ,  are  computed  without  taking  into 
1M  2M 

account  the  correlation  between  x^(t)  and  x2(t)  completely 

erroneous  results  are  obtained,  as  shown  m  Figure  1.  It 
should  be  noted  that  transfer  functions  are  complex  quantities 
in  general.  They  are  denoted  here  by  gain  and  phase  variables 

H ( f )  =  | H ( f )  j  ej  $  (f) 


-  4  - 


(7) 


Measured  System  Response  Functions 


where  j  H(f)  j  is  the  gain  factor  and  $  (f)  is  the  phase  factor. 

The  ordinary  coherence  functions  are  plotted  in  Figure 
2.  If  used  by  themselves  to  infer  the  nature  of  the  system  being 
investigated,  one  would  tend  to  conclude  that  high  measurement 
noise  (i.e.  system  noise)  was  present  or  that  nonlinear  effects 
were  present.  The  true  nature  of  the  system  is  exhibited  when 
the  partial  coherencies  are  computed,  as  shown  in  Figure  3.  The 
fact  that  all  three  functions  are  equal  to  one  over  the  frequency 
range  indicates  the  true  linear  relationships  which  exist  between 
X1  (t)  ,  x2(t)  '  and  * 

The  demonstration  cases  investigated  through  the  aid  of 
computer  simulation  methods  illustrate  the  usefulness  of  partial 
coherency  techniques  in  analyzing  multiple  input-multiple  output 
systems.  The  fact  that  correct  values  of  transfer  functions  and 
coherence  functions  are  obtained  when  all  in^.  its  are  properly 
taken  into  account  has  been  numerically  illustrated. 

The  next  step  in  the  investigation  will  be  to  use  the 
computational  procedures  to  compute  transfer  and  coherence  functions 
completely  from  actual  data.  These  can  then  be  compared  with  the 
expected  results  discussed  in  the  report.  Questions  of  statistical 
accuracy  which  will  be  limited  by  sampling  errors  in  the  spectral 
computations  can  then  be  explored.  Another  area  for  future  study 
is  the  trade-of^s  between  statistical  accuracy,  number  of  inputs, 
and  computation  time. 

Finally,  a  general  program  must  be  developed  for  the 
multi-input  case  and  applied  to  the  measurement  and  interpreta¬ 
tion  of  seismic  noise  in  large  arrays. 

B  Analysis  of  Arrival  Times  Across  TFSO  Extended  Array 

TFSO  extended  array  travel  time  anomalies  discussed  in 
the  15  April  1965  Semi-Annual  Technical  Summary  Report  has  been 
extended  to  include  events  from  six  teleseismic  epicentral  regions 
northwest  of  TFSO.  The  travel  time  anomalies  were  measured  from 
films. 


Table  I  shows  the  mean  travel  time  anomalies  relative 
to  the  center  of  TFSO  for  the  eight  LRSM  sites  in  the  extended 
array.  Some  of  the  events  used  to  measure  the  Kurile  Island 
anomalies  were  also  used  for  the  Honshu  and  Kamchatka  regions, 
however,  both  of  these  regions  include  measurements  from  events 
not  included  in  the  Kurile  Island  group. 
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Ordinary  Coherence  Functions 


Figure  2 
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Travel  Time  Anomalies  at  TFSO  Extended  Array 


The  size  of  the  travel  time  anomalies  do  not  seem  to  be 
consistent  for  all  six  regions  for  any  LRSM  station.  The  difference 
between  maximum  and  minimum  travel  time  anomalies  at  each  station 
varies  from  .12  seconds  at  GE  AZ  to  .75  seconds  at  NL  AZ.  The 
standard  deviations  in  the  travel  time  anomalies  are  consistent 
and  approximately  .10  seconds.  Moreover,  the  standard  deviations 
are  approximately  independent  cf  the  size  of  the  travel  time 
anomaly.  The  anomalies  from  the  four  central  regions,  which  are 
closest  to  the  same  azimuth,  are  more  consistent  than  results  from 
all  six  but  as  yet  no  obvious  functional  pattern  describing  the 
variation  has  become  evident. 

The  effect  of  summing  traces  in  the  extended  array 
without  some  kind  cf  station  corrections  is  shown  in  Figures  4 
and  5. 

Figure  4  shows  the  effect  of  computing  a  correlogram 
for  this  event  by  aligning  the  array  on  velocity  only.  The  first 
trace  is  the  velocity  summation  of  Z67,  GE  AZ,  SN  AZ,  Z-47,  and 
Z-74.  The  second  trace  is  the  (velocity)  summation  of  SG  AZ, 

JR  AZ,  LG  AZ,  WO  AZ,  and  NL  AZ. 

Figure  5  shows  the  correlogram  outputs  with  this  event 
aligned  by  velocity  plus  the  travel  time  anomalies  appropriate 
for  the  Kurile  Island  earthquakes  shown  above. 

A  method  has  been  developed  to  align  events  by  cross 
correlations.  The  array  is  aligned  approximately  by  velocity. 

Then  the  first  20  seconds  of  P-wave  for  each  trace  is  cross' 
correlated  with  the  similar  20  seconds  from  the  Z-47  trace.  The 
shift  of  the  peak  of  the  cross  correlation  function  from  the 
origin  indicates  the  required  shift  of  that  trace  relative  to  Z-47. 

The  results  using  this  method  ure  comparable,  if  not  better  than 
those  shown  in  Figure  5,  and  the  method  can  be  used  to  automate 
the  determination  of  station  corrections. 

C.  Analysis  of  Seismic  Data  From  WAGTAIL 

WAGTAIL  was  detonated  at  the  Nevada  Test  Site  on  3  March  1965. 
Thirty-four  stations  recorded  short-period  signals  and  20  stations 
recorded  long-period  signals.  It  is  possible  that  more  stations  in 
-he  U.  S.  would  have  recorded  surface  waves  if  it  were  not  for  the 
ins^rference  caused  by  an  unidentified  earthquake.  Thirteen  stations 
showed  compressional  first  motion  as  defined  by  the  First  Motion 
Criteria  (TWG  II).  The  average  magnitude  was  5.33  and  the  most 
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distant  station  recording  WAGTAIL  was  at  Oslo,  Norway,  a  distance 
of  about  8,130  kilometers.  Table  2  summarizes  the  measurements 
made  of  the  principal  phases  from  the  WAGTAIL  event.  Included 
are  the  Pn  and  P  arrival  times,  the  maximum  amplitudes  (A/T)  of 
Pn  or  P  and  Pg  motion  as  seen  on  the  shore-period  vertical 
instruments,  and  the  maximum  amplitudes  (A/T)  of  the  Lg  phases 
as  measured  on  the  short-period  horizontal  tangential  component. 
Long-period  Love  and  Rayleigh  wave  motion  are  also  tabulated  in 
(A/T)  form. 

D.  Computation  of  Depth  of  Focus  Using  pP  Arrival  Times 

The  epicenter  location  program,  LOCATE,  has  been  modified 
to  determine  depth  of  focus  using  pP  time  intervals.  After  having 
read  in  pP  arrival  times  for  an  event  in  the  same  manner  as  P, 
the  depth  is  determined  from  pP  by  first  computing  the  distance  from 
tie  approximate  epicenter  to  all  stations  reporting  pP  times;  the 
table  of  pP-P  time  interval  vs.  distance  is  entered  to  find  the 
estimated  depth  of  focus  for  each  station  reporting  pP  time;  the 
mean  and  standard  deviation  of  depth  is  computed  from  station 
depth  estimates;  finally,  epicenter  and  origin  time  are  computed 
with  depth  constrained  to  the  computed  value  as  determined  by 
the  prior  computations.  Upon  convergence  of  the  program,  the 
confidence  regions  are  computed  as  before  for  latitude,  longitude 
and  origin  time.  In  the  arrival  time  listing,  the  pP  and  pP-P 
intervals  are  printed  along  with  the  other  information. 

E.  Analysis  of  Short-Period  Noise  and  Arrays  at  Three 

English  Seismological  Observatories 

This  study  was  made  jointly  by  the  VELA  Seismological 
Center  (VSC)  and  the  Sei '  ' ' c  Data  Laboratory  and  resulted  in 
AFTAC  Technical  Report  VU-o5-2. 

Programs  were  written  and  the  following  data  computations 

made: 


1.  Five  150  second  multichannel  data  samples  were 
digitized  at  20  samples/sec  to  produce  3000  data 
points  from  each  channel.  The  mean  and  linear 
trend  were  removed. 

2.  High  resolution  (300  lags)  and  low  resolution  (100 
lags)  power  spectra  were  computed  for  one  channel 
from  each  multichannel  sample. 
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3.  After  filtering  and  decimating,  geometrically 
non-redundant  cross  correlations  were  computed. 

These  were  transformed  in  frequency  and  wave 
number  to  produce  wave-number-frequency  (w-n-f) 
power  spectra. 

4.  Plots  of  all  power  spectra  and  (w-n-f)  power 
spectra  were  produced,  together  with  plots  of 
the  infinite  velocity  summation  response  and 
the  transform  of  the  spatial  window  for  the 
(w-n-f)  spectra. 

F.  Seismic  Data  Compilation  and  Retrieval 

The  purpose  of  this  study  was  to  prepare,  from  several 
programs  previously  written,  a  program  to  produce  compilation  of 
the  seismic  activity  of  the  earth  by  regionalization. 

The  program  has  been  written  and  is  operational;  seismic 
data  for  the  years  1960  through  1964  has  been  compiled  for  one-by- 
one,  two-by-two,  and  five-by-five  degree  blocks. 

The  program  has  the  capability  to  compute  a  function  of 
magnitude,  F(M) ,  for  all  events  in  any  given  year  or  years,  the 
function  being  summed  over  blocks  K  by  L  degrees  in  size,  where 
K  is  a  divisor  of  180  and  L  is  a  divisor  of  90.  The  function 
K(M)  is  at  present  simply  a  count  of  events,  but  it  can  be 
replaced  by,  e.g. ,  number  of  events  per  unit  area,  number  events 
smaller  or  larger  than  a  given  magnitude,  etc.  Also  the  program 
can  be  run,  using  functions  ether  than  a  simple  count. 

III.  WORK  IN  PROGRESS 

A.  Teleseimic  Signal  and  Noise  Correlations  at  the  TFSO 

Extended  Array 

The  objective  of  this  study  is  to  determine  the  largest 
dimensions  of  a  seismic  array  for  which  the  ordinary  array  processing 
of  aligning  signals  and  mixing  (i.e.,  velocity  filtering)  still 
reinforces  the  signals  and  cancels  the  noise.  Velocity  filtering 
tends  to  enhance  signals  of  the  same  shape  and  tends  to  smear, 
average,  or  cancel  signals  of  different  shapes.  The  correlation 
coefficient,  which  is  the  maximum  value  of  the  normalized  cross 
correlation  function,  provides  a  quantitative  measure  of  the 
similarity  in  waveshape  between  two  seismic  signals.  Plotting 
the  correlation  coefficients  vs.  separation  between  seismometers 
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will  provide  a  measure  of  correlation  vs.  distance  over  the  array, 
and  filtering  the  seismograms  through  narrow  bandpass  filters 
-»  correlating  will  provide  a  measure  of  correlation  vs. 
freq  cy. 


For  this  study,  the  signals  used  were  P-waves  from  10 
teleseismic  earthquakes  recorded  over  the  TFSO  extended  array. 

Noise  samples  u3ed  were  from  these  same  events  just  prior  to  the 
onset  of  the  signals.  The  events  themselves  were  chosen  for  the 
largest  signal-to-noise  ratio  possible  which  avoided  clipping  on 
the  magnetic  tapes. 

B .  Correlogram  Analysis  from  Linear  Arrays 

Corre'ogram  analysis  is  a  continuing  study,  testing  the 
British  array  processing  technique.  Earthquake  data,  recorded  at 
5  to  10  array  stations,  for  the  months  of  May,  June,  July,  August, 
and  December  1963,  January,  February,  and  March  1964,  have  been 
processed  and  data  for  April  1964  are  being  processed. 

C .  Dispersion  Analysis  of  Surface  Waves  from  Deep  Well  Data 

In  addition  to  the  10  deep  wells  previously  studied  and 
reported,  a  study  of  the  seismic  noise  recorded  in  one  deep  well 
was  made.  The  analysis  included  est:'  lates  of  power  spectral  density, 
phase  angle,  and  coherency,  and  the  theoretical  Rayleigh  waves  in 
the  structure. 

P.  Rectilinear  Motion  Detection  (REMODE) 

The  detection  and  measurement  of  linaarly  polarized 
particle  motion  is  underway.  Various  methods  for  detecting  time 
varying  polarized  waves  are  described  in  the  literature.  These 
along  with  our  own  ideas  are  being  evaluated  in  order  to  improve 
the  detection  process.  The  following  questions  are  being  asked: 

1.  What  S/N  gain  is  obtained  for  strong  signals? 

2.  What  is  the  input  S/N  threshold  for  detecting  weak 
signals? 

3.  How  badly  is  the  signal  being  distorted  by  the 
process? 

The  data  for  answering  these  questions  experimentally 
consists  of  the  following: 

1.  A  strong  signal  with  virtually  no  noise  is  recorded. 
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2. 


The  signal  Z  and  R  component  are  added  to  noise 
recorded  at  the  same  station. 

3.  The  noise  level  is  increased  in  steps  until  the 
signal  is  not  detectable  by  the  process. 

The  following  measurements  are  manually  made  using  the 
above  experimental  data 

1.  Peak  to  peak  amplitude  of  the  signal  plus  noise 
(including  the  signal  alone  and  noise  alone)  is 
taken  in  the  time  window  of  the  signal's  maximum 
peak  to  peak  amplitude . 

2.  These  measurements  are  used  to  compute  the  S/N 
ratio  of  the  unprocessed  and  processed  record. 

The  S/N  of  unprocessed  records  is  plotted  vs.  the 
S/N  of  processed  records.  Both  of  these  quantities 
are  based  on  an  average  computed  from  measurements 

of  the  signal  added  to  twenty  different  noise  samples 
at  the  sam  !  preset  noise  level.  The  noise  level  is 
changed  approximately  12  times,  sufficiently  to  draw 
precise  experimental  curves . 

The  results  to  date  have  been  used  to  evaluate  the  method 
of  Shimshoni  and  Smith,  using  the  scalar  correlation  coefficient  as 
a  gain  control.  Approximately  6  to  12  db  of  gain  was  obtained  by 
the  process.  The  detection  threshold  was  for  signal  approximately 
0  db  above  the  noise. 

The  same  method  was  evaluated  with  the  modification  of 
rotating  the  Z  and  R  component  to  Z'  and  R'  determined  so  that  the 
components  are  nearly  equal  (as  contrasted  with  an  initial  orbit, 
R/Z,  of  .25;  R'/Z'  =  1.0).  This  produced  an  increase  in  the  gain 
of  6  to  12  db,  so  that  18  db  was  obtained  by  the  process  for  strong 
signals.  Furthermore  the  detection  threshold  was  pushed  back  to 
signals  2  db  a^ove  the  noise. 

So  far  all  the  results  have  been  obtained  on  the  analog 
computer.  This  work  is  continuing.  The  testing  procedure  is  being 
extended  co  the  digital  computer  ard  the  digital  program  used  to 
process  the  Salmon  data  plus  any  refinements  to  it  (such  as  the 
rotation  method  previously  described)  will  be  similarly  evaluated. 

E .  Automation  of  Data  Processing  and  Analysis 

Preliminary  processing  of  time  series  data  has  become 
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relatively  standardized  at  SDL.  In  particular,  operations  of 
filtering,  spectral  analysis,  auto-  and  cross-correlation, 
measurements  of  coherency,  instrument  response  and  elementary 
summing  and  multiple  ition  of  time  series  are  used  repeatedly, 
often  many  of  these  operations  being  carried  out  on  the  same 
set  of  data.  In  addition, many  standard  special  purpose 
analysis  programs  are  now  in  use  and  the  number  of  such  programs 
is  expanding  at  a  rapid  rate. 

It  is  advantageous  then  to  consider  combination  and 
automation  of  these  processing  and  analysis  methods,  especially 
those  that  have  already  become  standardized  or  may  reasonably  be 
s' andardized  in  their  function.  In  particular,  generalized  multi¬ 
purpose  chain  programs  consisting  of  several  specialized  processing 
and  analysis  programs  can  serve  to  combine  these  often  lengthly 
procedures  and  to  thereby  increase  the  efficiency  of  data  processing 
to  the  point  where  more  complicated  operations  become  feasible. 

A  system  program  FLAP  (FORTRAN  Controlled  Linkage  of 
Associated  Programs)  has  successfully  been  designed.  Basically, 
the  system  functions  rs  a  subroutine  control  system  which  allows 
a  FORTRAN  program  (executive  program)  to  call  and  execute  programs 
in  the  same  manner  as  subroutines.  The  system  fulfills  all  the 
requirements  outlined  above  with  added  features  of  flexibility 
inherent  in  a  purely  FORTRAN  system.  In  particular,  digital 
operations  (in  the  FORTRAN  language)  may  be  programmed  between 
calls  to  the  individual  programs.  This  results  in  optimal 
interfacing  between  the  various  special  programs. 

A  few  chain  systems  involving  seismic  wave  analysis 
and  processing  are  now  in  use  and  many  others  are  in  various 
stages  of  completion.  Such  standardized  package  programs 
(chains)  provide  well  defined,  processed  data  with  a  minimum 
of  time  and  effort. 

F.  Recursive Numerical  Filters 


A  set  of  programs  which  can  be  used  to  model  analog 
filters  on  the  digital  computer  are  now  operational.  Certain 
problems  still  exist  in  the  use  of  these  programs,  but  sufficient 
checks  are  built  in  so  that  the  inherent  errors  should  not  be 
a  serious  limitation  in  practice. 

Seismic  systems  have  been  modeled.  The  major  emphasis, 
however,  has  been  on  the  design  of  low-  and  band-pass  filters. 
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For  these,  the  program  limitations  have  been  well  bounded. 

Stable  low-pass  filters  have  been  designed  with  roughly 
the  following  limits:  low-pass  band  width  of  0.10  cps.  with 
24  db  slopes,  0.30  cps.  with  36  db  slopes,  1.0  cps  with  60  db 
siopes,  and  8.0  cps.  with  120  db  slopes. 

The  band-pass  filters  are  designed  by  first  generating 
the  coefficients  of  the  low-pass  filter  which  is  half  as  wide 
as  the  requested  band-pass  filter,  and  then,  using  the  coefficients, 
a  new  set  of  coefficients  is  computed  by  shifting  the  low-pass 
center  frequency  from  0.00  cps  to  the  requested  center  frequency. 

To  date,  no  stable  filter  has  been  nu'de  unstable  by  shifting. 

At  present,  it  is  felt  that  no  further  research  type 
work  should  be  done  on  the  program.  However,  as  the  program  now 
exists  in  independent  parts,  the  present  plan  is  to  put  the 
parts  together  and  introduce  input  flexibility  so  that  the 
program  irny  be  presented  as  a  package  to  the  general  user. 

G.  Modification  of  LOCATE  Program  for  Event  Epicenter 

Determination  using  Reported  P  Arrival  Times  from 

3  or  4  Stations. 

The  event  location  program,  LOCATE,  has  been  revised 
to  accept  data  from  a  minimum  of  three  stations.  With  this  small 
number  of  stations,  the  standard  errors  ar d  confidence  regions  can 
not  be  computed  as  they  are  normally  done  when  data  from  more 
than  five  stations  are  used.  Moreover,  to  operate  on  data  from 
three  stations,  two  of  the  four  normal  variables  are  constrained 
to  the  input  values;  these  are  the  origin  time  and  depth  of  focus, 
aid  are  taken  from  the  USC  &  GS  listing  which  are  assumed  to  be 
correct  foi  each  event. 

In  an  operational  system,  the  origin  time  and  depth 
would  not  be  known;  in  which  case,  the  depth  would  be  assumed  to 
be  zero  or  some  other  convenient  value. 

A  preliminary  study  has  been  made  to  test  the  precision 
of  epicenters  based  on  3  and  4  station  data. 
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Only  relatively  large  magnitude  events  were  used,  with 
data  from  a  very  few  falling  in  the  low  magnitude  four  range. 

The  results  ,-,ere  studied  in  terms  of  the  computed  azimuth  and 
distance  of  each  event  from  the  epicenter  determined,,  using  data 
from  a  la  : je  number  of  well  distributed  stations,  as  reported  by 
the  USC  &  GSo 

The  attached  Table  3  lists  the  epicenter  shifts  in 
kilometers  and  azimuth  along  with  the  standard  deviations  and  the 
measurement  of  dispersion. 

We  can  conclude  from  this  limited  data  that  an  epicenter 
based  on  data  from  three  or  four  stations  is  valid  within  approximately 
200  kilometers.  Although  this  precision  is  obviously  not  good 
enough  to  help  us  make  a  practical  decision  regarding  on-site 
inspection,  it  is  good  enough  as  a  first  approximation  on  which 
we  may  base  corrections  and  techniques  for  precise  beam  steering 
and  a  precise  epicenter  determination  based  on  data  from  many 
stations . 


The  tripartite  array  TFO,  UBO,  and  BMO  is  roughly  within 
the  same  geologic  provice  while  WKO  in  Oklahoma  clearly  causes 
the  array  TFO,  UBO,  and  BMO  to  overlap  two  distinct  provinces. 

This  fact  may  explain  in  part  the  apparent  slight  improvement  of 
epicenter  computations  based  on  the  first  array  over  those  based 
on  the  second  array.  The  use  of  all  four  stations  appears  to 
improve  the  precision  even  more,  but  probably  not  enough  to  improve 
appreciably  the  choice  of  arrival-time  corrections  needed  for  beam 
steering  and  precise  epicenter  determination.  The  principal 
advantage  of  four  stations  is  provide  redundancy  in  the  event  that 
data  from  one  of  the  stations  is  unusable. 

Further  work  on  this  problem  should  include  similar 
studies  of  events  from  other  regions  of  the  world,  the  effect  of 
magnitude  or  signal-to-noise  ratio  on  epicenter  determination,  and 
the  effect  of  making  station  travel  time  corrections. 

H .  Large  Aperture  Seismic  Array  (LASA)  Data  Processing 

Studies  are  being  made  on  the  additional  electronic 
equipment  and  the  new  data  processing  requirements  which  the  SDL 
will  need  in  order  to  process  the  LASA  data. 

Equipment  ordered  includes  a  604  magnetic  tape  drive  and 
controller  which  will  read  LASA  data  generated  at  800  bpi,  a 


13 


TABLE  OF  COMPUTED  ERRORS 
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200  million  bit  disk  file  which  will  afford  rapid  access  to  the 
LASA  data.  Still  under  consideration  are  several  types  of  digital- 
to-analog  su^ystems  for  data  processing. 

Two  digital  computer  programs  have  been  written.  One 
program,  written  and  checked  out  for  the  160-A,  will  read  LASA 
data  and  transmit  it  from  the  160-A  to  the  1604  via  the  satellite 
hardware.  The  other  program,  written  for  array  processing,  receives 
the  data  from  the  160-A,  phases  the  data  for  a  specific  angle, 
velocity,  starting  time,  and  time  interval,  and  outputs  a  tape 
containing  the  phased  sum  of  each  of  the  21  subarrays. 

A  maximum  of  20  data  points  can  be  lagged  across  each 
subarray.  Multiple  input  and  output  tapes  are  allowed  and  any 
amount  of  input  data  can  be  processed.  The  array  configuration 
can  be  varied.  It  is  estimated  that  processing  will  take  4.5 
times  real  time. 

IV.  SUPPORT  AND  SERVICE  TASKS 

As  part  of  the  contract  work-statement,  the  SDL  provided  one 
or  more  of  the  following  support  and  service  functions  for  AFTAC 
and  other  VELA  participants: 

-  copies  of  16  and  35  mm  film 

-  copies  of  existing  composite  analog  tapes 

-  composite  analog  tapes  of  special  events 

-  use  of  1604  computer  for  checking  out  new  programs 
or  running  production  programs 

-  copies  of  digital  programs 

-  digitized  data  in  standard  formats  or  special  formats 
for  use  on  computers  other  than  the  1604 

-  running  SDL  production  programs,  such  as  power  spectral 
density,  and  array  processing  on  specified  data 

-  digital  x  -  y  plots  of  power  spectra  or  digitized  data 

-  signal  reproduction  booklets 

In  addition,  visiting  scientists  utilized  SDL  facilities  to 
study  data  and  exchange  information  with  SDL  personnel. 

A.  Data  Library 

The  Data  Library  contains  3600  digitized  seismograms,  106 
digital  computer  programs,  and  212  composite  analog  magnetic  tapes, 
all  available  for  use  by  the  VELA-UNIFORM  program. 
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period: 


The  following  additions  were  made  during  this  report 


1.  Digital  Seismograms  -  554,  including: 

-  deep  well  noise  from  EK  NV,  OR  FL,  and  AP  OK 

-  noise  samples  on  a  new  slow-speed  tape 
transport 

-  long  and  short  period  data  from  BILBY,  SHOAL, 
and  the  FALLON  earthquake 

-  noise  data  from  one  of  the  LASA  subarrays 

-  many  earthquakes  recorded  at  the  VELA 
observatories  and  the  TFO  extended  array 

2.  Digital  Programs  -  4,  including: 

-  POLES1A2  -  This  program  generates  the  complex 
plane  poles  and  normalizing  constant  for  use 

in  characterizing  the  denominator  and  numerator, 
respectively,  of  a  LaPlace  transform,  /it  present, 
the  subroutine  creates  only  Butterworth  or 
Chebyschev  poles.  However,  it  is  constructed 
to  be  easily  expanded  to  include  other  special 
functions. 

-  RESPOND  -  For  given  wave  number,  this  program 
computes  the  response  of  an  array  of  seismometers 
along  a  specified  azimuth. 

-  DCONVOL1  -  To  remove  the  response  in  a  seismogram 
due  to  the  inherent  filter  characteristics  of 
the  seismometer,  filter,  velocity  transdxicsr, 

and  galvanometer.  The  routine  removes  the 
response  of  one  instrument  at  a  tine,  using  the 
filter  characteristics  of  that  instrument. 
Successive  passes  through  the  routine  with 
different  filter  characteristics  will  remove 
all  the  undesired  responses. 

-  LAPINV  -  Given  a  function  of  f'  e  complex 
variable  s,  and  an  array  of  the  independent 
variable  time,  this  program  generates  the 
time  function  for  each  of  the  time  values 
indicated. 
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programs ^ 


In  addition,  modifications  were  made  to  the  following 


a.  The  LOCATE  program  has  been  modified  to  provide  the 
option  of  computing  depth  ot  focus  from  pP  arrival 
times,  and  also  to  locate  events  using  3  or  more 
stations,  instead  of  a  minimum  of  6= 

b.  The  COLLATE  program  has  been  modified  to  run  from 
the  normal  epicenter  tapes  instead  of  a  special 
merged  tape  as  it  required  before.  This  modification 
increases  the  efficiency  of  the  program  and  reduces 
the  running  time  by  a  factor  of  10. 

c.  The  DEPTHMAG  program  has  the  following  added  features 

(1)  ability  to  retrieve  from  any  one  of  four 
differently  formated  tapes. 

(2)  ability  to  save  all  card  images  retrieved. 

(3)  ability  to  specify  retrieval  criteria  relative 
to  recording  stations,  data  source  codes,  and 
distance  and  azimuth. 

3.  Analog  Composite  Tapes  -  20,  including: 

a.  Made  by  Seismic  Data  Laboratory 

-WAGTAIL 

-Special  composites  made  for  AFTAC,  the  Air 
Force  Cambridge  Research  Laboratory,  Texas 
Instruments,  Inc,,  and  General  Atronics 
Corporation 

b.  Made  by  the  Geotechnical  Corporation 

-TURF 

-STURGEON 

-OCONTO 

-HOOK 

-TORNILLO 

4 .  Other  Data  Available  to  VELA  Participants 

The  following  data  was  added  to  our  magnetic  tape  storage: 
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-  U.  S.  Coast  and  Geodetic  Survey  epicenters 
from  March  1965  to  May  1965 

-  Earthquake  Bulletin  Data  from  the  LRSM  stations 
(November  1964  to  February  1965) 

-  Earthquake  Bulletin  Data  from  the  VELA 
observatories  (December  1964) 

B.  Data  Compression 

This  is  a  continuing  routine  operation,  and  production  is 
maintained  at  the  level  needed  to  meet  the  requirements  of  the  field 
operations  (LRSM  teams  and  U.  S.  observatories)  and  the  Seismic  Data 
Laboratory. 

C.  Equipment  Modifications 


No  modifications,  other  than  those  discussed  under  III-H 
(LASA  planning) ,  were  made  during  this  report  period. 

D.  Automated  Bulletin  Process 


During  this  report  period,  the  ABP  digital  program  was 
run  for  the  January,  February  and  March  1965  earthquakes  bulletins 
issued  by  the  Geotechnical  Corporation. 

The  program  is  being  continuously  modified  to  associate 
and  identify  a  greater  percentage  of  phase  arrivals. 


ORGANIZATIONS  RECEIVING  SDL  DATA  SERVICES 


APRIL  1965  -  JUNE  1965 


AFCRL 

Boston  College 
Engineering  Physics 
General  Atronics  Corp. 
Geotechnical  Corp. 

Institute  for  Defense  Analysis 
International  Seismological 
Research  Centre ( Scot land) 


Lamont  Geophysical  Observatory 
Lincoln  Laboratories 
Penn  State  University 
Princeton  University 
St.  Louis  University 
Texas  Instruments,  Inc. 

U.S.  Coast  and  Geodetic  Survey 
Vitro  Corp. 
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